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large super-Hubble modes. 
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I. INTRODUCTION 

The unprecedented isotropy of the cosmic microwave background radiation revealed by 
modern CMB experiments is - justifiably enough - usually taken as an indication that 
we inhabit an isotropic universe. Together with the apparent (but far more debatable) 
homogeneity of matter distribution at scales > 100/i^^Mpc, this observation suggests a 
generalization of the Copernican principle: our universe is, on average, homogeneous and 
isotropic. At present, the Copernican principle remains a hallmark of cosmology, and its 
modification would profoundly alter the way we understand the physical universe. As a 
matter of fact, challenges to this scenario have largely boosted research in cosmology in the 
past, and will certainly do so in the foreseeable future. 

In the last decade, new possibilities for a different approach to the cosmological principle 
have emerged from two different types of observations. On the one hand, the detection of 
the accelerated expansion of the universe [1-3] has motivated dark-energy-free cosmological 
models in which our (allegedly) privileged position in the universe would be perceived as 
a global cosmic acceleration [4, 5]. On the other hand, discovered large-angles Cosmic 
Microwave Background (CMB) anomalies are being interpreted as a manifestation of a 
preferred direction in the universe (see Ref. [6] for a recent review). Regardless of which 
description might turn out to be the best, these investigations are certainly well founded 
from the observational point of view - not to mention that they constitute a necessary step 
to be faced in the so-called era of precision cosmology. 

From a theoretical perspective, it is legitimate to ask when one can establish a one-to-one 
relation between our observations and the symmetries of the universe we live in. Although 
a homogeneous and isotropic universe would certainly lead to the observed CMB, is the 
converse true? Actually, Ehlers, Geren and Sachs (EGS) proved that the converse holds in 
an expanding universe provided that the isotropic CMB is measured by fundamental ob- 
servers in a dust universe [7]. However, there are important counterexamples which allow 
the engineering of a richer class of models without jeopardizing well-established cosmolog- 
ical observations [8-10]. The theorem demonstrated by EGS admits generalizations to a 
matter content given by a generic perfect fiuid [8], and for an almost isotropic CMB [11], 
case for which the associated geometry is only approximately homogeneous and isotropic. 
But, how can we phenomenologically conceive an anisotropic model of the universe in which 



basic observables, like for instance the redshift, remain isotropic, at least at dominant or- 
der? One example was proposed by Minioso and Crawford in 1993 [12]. They showed that 
some anisotropic spacetimes, namely Bianchi type-III (Bill) and Kantowski-Sachs (KS) 
spacetimes, admit a shear-free expansion in which the scale factor has the same dynam- 
ics as a spatially curved Friedmann-Robertson- Walker (FRW) universe, provided that the 
anisotropic stress-tensor of the cosmological budget is in direct proportion to the electric 
part of the Weyl tensor. Since the dynamics is governed by a single scale factor, the met- 
ric can be brought to a conformally static form. It then follows by another EGS result 
that the radiation will be isotropic, in accordance with basic observations [7] (see also [10]). 
Further and independent developments have shown that this scenario can also be realized 
by endowing the universe with other forms of matter fields. In Ref. [13], for example, a 
shear-free model with anisotropic curvature was constructed by means of a two-form field. 
Another interesting possibility was proposed in Ref. [14] (and further discussed in Ref. [15]), 
where a massless and anisotropic scalar field was used to counterbalance the curvature of a 
Bill universe, producing an anisotropic model with a FRW-like background dynamics. This 
latter example is particularly interesting, given the role that scalar fields play in modeling 
early-universe physics. 

Evidently, this alternate scenario may become appealing and competitive only if one is 
able to extract distinctive cosmological signatures from it. A possible observable at the 
level of the background geometry is the luminosity distance function: because the geome- 
try traversed by photons along our line of sight is characterized by an anisotropic (spatial) 
curvature, the luminosity distance function will not have a direct correspondence to the 
standard one used in FRW models. One can then envisage a setting in which supernova sur- 
veys can be used to detect anisotropic modulations of the luminosity function [13], therefore 
rendering the proposed model falsifiable. However, if one wants to confront this theoretical 
scenario against cosmological data, then a much broader class of observables can be attained 
by considering perturbation theory, rather than focusing on the background. In particular, 
the whole framework of the inflationary era and its connection to structure formation pro- 
vides a vast arena within which one can discuss the physical consequences of anisotropic 
models with isotropic expansion and put them to the test. 

The theory of linear perturbations in non-FRW cosmological backgrounds has been an 
active field of research by its own, where many different models have been studied. The 



case of isotropic and inhomogeneous Lemaitre-Tolman-Bondi spacetinies, for example, has 
been recently worked out in Ref. [16], while the theory of perturbations in homogeneous and 
globally anisotropic Bianchi I spacetimes has been presented in Refs. [17-19]. Quite inde- 
pendently, the dynamics of linear perturbations in spatially curved anisotropic backgrounds 
has been investigated in Ref. [20] , where the Bill case has been considered, and also in Ref. 
[21], where both the Bill and KS cases have been analyzed. 

This paper presents the formalism necessary to extract the signatures that shear-free 
anisotropic models of the Bill family would generically produce in linear cosmological per- 
turbations, specialized to situations where the universe is dominated by a slowly-varying 
scalar field, as it typically happens during inflation. With this aim, as well as to make more 
precise statements such as "up to conditions on the behavior at infinity", which are usually 
employed in a rather vague sense, we introduce an explicit set of complete eigenfunctions 
in terms of which perturbations can be expanded. This is a crucial step in perturbation 
theory that is sometimes overlooked in the literature, and unavoidable if we want to apply 
perturbation theory to specific cosmological models. To the best of our knowledge, these 
eigenfunctions have not been considered elsewhere. 

As we will show, during a generic slow-roll inflationary phase, the dynamics of the per- 
turbations can be described by effective degrees of freedom which have a remarkable re- 
semblance to the standard Mukhanov-Sasaki (MS) variables [22, 23]. However, since the 
eigenmodes of the Laplacian on the spatial sections are not given by standard plane-waves, 
we might expect non-trivial features in the morphology of both the temperature and matter 
fields of the universe. A particular signature which arises immediately from the spectrum 
of the Laplacian is the obstruction to have arbitrarily large super-Hubble modes. Because 
the curvature naturally introduces a spatial scale, the largest possible wave modes will be 
comparable to this scale. We also comment on the possibility of detecting specific features 
of shear-free anisotropy in the spectrum of gravitational waves. 

This paper is structured as follows. We begin by reviewing the fundamental ideas behind 
anisotropic shear-free spacetimes in Sec. II, where we also introduce basic notations and 
definitions. In Sec. Ill we discuss the appropriate mode decomposition of the perturbations 
and use it to construct gauge-invariant variables. We close this section by determining a 
convenient basis of eigenfunctions of the Laplace operator. The dynamical equations of the 
cosmological perturbations are presented in Sec. IV, where we also construct the anisotropic 



equivalent of the MS variables. We conclude in Sec. V. In this work we will adopt units such 
that c = 1 = SttG, and spacetime signature (— , +, +, +). We follow the standard convention 
that Greek and Latin indices label spacetime and spatial coordinates, respectively. The 
spacetinies we will consider have a residual symmetry along one direction. We choose this 
preferred direction to lie along the 2;-axis and let the indices {a, b, c} run from 1 to 2, so as 
to span the xy-plane. 



II. SHEAR-FREE BACKGROUND DYNAMICS 

The simplest example of anisotropic universe is the one described by a Bianchi I spacetime. 
Since in this family of models the shear is entirely due to the anisotropic expansion of the 
universe, the spacetime cannot be both anisotropic and shear-free. This behavior is of 
course expected, since Bianchi I models contain the fiat FRW case in the shear- free limit. 
The possibility of having anisotropic shear-free models with isotropic expansion relies on 
the fact that some anisotropic spacetimes do not possess a FRW limit. In these cases, the 
anisotropy of the space does not have to be encoded in the shear, but can instead be caused 
by imperfections of other kind in the fluid filling the universe. 

An example of such spacetime is found in the Bill type. The Bill spacetime in question 
is spatially curved and possesses residual isotropy with respect to one direction. We choose 
the preferred axis of symmetry to point along the z direction. Spatial sections are described 
by the direct product of the real line with a pseudo-sphere, M x H^, thus providing a space 
with negative constant curvature. In this work we will adopt a coordinate system in which 
the shear-free Bill metric can be written as: 

ds^ = air^f {-drf + Tabdx^'dx' + d^^) ^ (1) 

where a{rj) is the isotropic scale factor, rj = J dt/a is the conformal time, and jab is the 
standard metric of EP, which in a convenient coordinate system x"" = {x,y} reads: 

7abdxMx'' = dx^ + e^'^dy^ . (2) 

Following [14], we will consider a matter sector composed by the standard infiaton field ip 



plus a massless scalar field ^: 

K) u = d'vd.^ - ^t (^^V^A^ + V{v)^ , (3) 

^W . = d^^d,^ - K {]^dHdxA ■ (4) 

The inflaton field being homogeneous, ip = ip{ri), satisfies the standard Klein-Gordon equa- 
tion: 

ip" + 2n^' + aV^ = . (5) 

The prime denotes the derivative with respect to the conformal time. The massless field, on 
the other hand, obeys 

^" + 2H'd' - V^^ = . (6) 

Here, V^ = D'^Da + <9| and Da is the covariant derivative compatible with 7^^. It is easy 
to check that the choice d = Cz, where C is an arbitrary constant, not only solves Eq. (6), 
but also produces a diagonal energy-momentum tensor: 

pw = pU = -Pw = ^ ' 

where pi^^ and p[L represent the pressures perpendicular and parallel to the xy-plane, re- 
spectively. With these definitions, we can proceed and obtain Einstein equations: 

^2 

(7) 

(8) 

(9) 

It is then clear that if we choose C^ = 1,^ Eqs. (8) and (9) turn out to coincide, while 
the whole system of equations becomes the standard Friedmann equations with curvature 
K= -1/2: 

H^ = ^^ + aV(v.))-/., (10) 

H' = -li^" - a'V{^)) . (11) 



37/2 = 


a^PH+ 2+1' 


2 + 2n' = 


q2 


2 + 2n' = 


(j2 



^ The ambiguity in the sign of C comes from the freedom in orienting the z-axis, which can be chosen at 

will. 
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This is indeed remarkable, since the dynamics is essentially that of an isotropic FRW 
universe with a perfect fluid source. On the other hand, it is worth noticing that our 
discussion surpasses the model obtained with our specific choice of solution for the scalar 
field; indeed, the same results can be obtained by endowing the universe with matter fields 
sustained by other considerations, such as a two-form field [13] or an imperfect fluid with 
stress [12]. In this work we will adopt a rather pragmatic route and consider the solution d = 
Cz as a way to unveil the generic signatures of inflationary perturbations in an anisotropic 
spacetime, while experiencing isotropic dynamics at the background level. As we will see, 
the main signatures of an inflationary era are rather insensitive to the perturbations of this 
field. 

In order to find reduced dynamical variables for the perturbations, it will prove convenient 
to work with Eqs. (10) and (11) combined in the following form: 

^-1 + 2^-^" (12) 

Finally, let us also introduce here the standard slow-roll parameters, which will be needed 
to identify the effective degrees of freedom during inflation. They are defined by: 

The first order slow-roll approximation is defined by the requirements that e and 6 remain 
small during inflation ■^: 

e < 1 , (5 < 1 . (14) 

This is always true provided that the inflaton evolves under a considerably flat potential. In 
this regime, one can also show that: 

n^---^, (15) 

r]l — e 
where, we remind the reader, tj = J dt/a = J da/^aTi). 

III. LINEAR PERTURBATION THEORY 

We start by describing the general framework for the consideration of linear perturbations. 
The metric and its inverse are written as a background piece plus contributions treated as 
^ The constancy of e and S is automatically ensured because their time derivatives are of second order in 
these parameters [24]. 



perturbations: 

9^..^9,. + Sg,,, gf^'^ ^ gf^-^ + Sg^^\ (16) 

At linear order in perturbation theory, this iniphes that 



To manipulate indices in perturbed quantities, one should apply (16) and then discard second 
order terms in the perturbations. 

The matter sector is composed by the fields {}, (p, and their perturbations, which can be 
written as: 

i^^^ + 6^, (17) 

^ -^ (p + 6ip. (18) 

In the following subsections we will discuss the general parametrization of 6g^i, and the 
construction of gauge-invariant variables. 

A. Mode decomposition 

The general construction of perturbation theory in cosmology is usually based on the 
identification of a set of independent modes representing the degrees of freedom of the 
perturbations. For spacetimes admitting a 1+3 foliation, these modes are typically iden- 
tified with irreducible pieces of the symmetry group of the corresponding 3-dimensional 
spatial sections. In the FRW case, these are maximally symmetric spaces, and a unique 
decomposition of symmetric rank-2 tensors in terms of Scalars, (transverse) Vectors and 
(transverse-traceless) Tensors is always possible (SVT decomposition, from now on) [25]. In 
the cases studied here, the constant-time hypersurfaces are not maximally symmetric spaces, 
so the standard SVT decomposition is not adequate. However, since the considered Bill 
spacetimes have a 2-dimensional maximally symmetric subspace, namely the pseudo-sphere, 
we can instead proceed with a 1+2+1 decomposition of the dynamical modes. 

Some important differences with the standard SVT decomposition are immediate: given 
that the largest subspace is of dimension 2, independent transverse and traceless symmetric 
tensor modes cannot exist on it, and the most general symmetric rank-2 tensor is decomposed 
as 

Tat = X^al> + DaD,Y + D^aZt), I?"^a = , (19) 

8 



where X and Y are scalars, Za is a transverse vector, and we recall that Da is the covariant 
derivative in H^. Likewise, any vector Va will be decomposed as 

Va = DaV + Va, I^'^K = , (20) 

where V^ is a scalar and Va is again transverse. Notice that the transversality condition 
eliminates one of the two variables of the transverse vectors, so that they are described by 
one single scalar function. Nonetheless we will maintain the standard terminology and refer 
to the splitting (19) and (20) as a Scalar- Vector (SV) decomposition. In the following, every 
barred vector will refer to a transverse vector. 

Based on the discussion above, we can parametrize the most general perturbation 5g^u 
as 

ds^ = a^ [-(1 + 2A)drf + 25„da;"dr/ + 2E„dx"d;z + 2Cdr]dz + /i^fedxMx^ + 2Fdz^] , (21) 

where 

Ba = DaB + Ba , (22) 

Ea = DaE + Ea , (23) 

hab = 2S^ab + 2DaDbU + 2D(,H) • (24) 

This decomposition is unique, and the introduced functions and vectors completely specify 
the 10 degrees of freedom in 6gfj_u. 

B. Gauge-Invariant Variables 

Let us specify the gauge transformations of the above free functions. Under an infinites- 
imal gauge transformation of the form 

xf' = x^'-C, (25) 

any perturbed quantity 6Q transforms as 

6Q ^ 6Q = 6Q + C^Q, (26) 

where Q denotes the background quantity to be perturbed. By definition, gauge-invariant 
variables are those for which 

/:^Q = 0. (27) 



According to the 1+2+1 splitting of the spacetime, the gauge vector ^^ will be parametrized 
by 

where T, L, and J are scalars and L^ is a transverse vector. 

1. Metric 

If the quantity Q is the spacetime metric, then the gauge transformation (26) becomes 

where V is the full covariant derivative compatible with the background metric (1) and, as 
usual, the parentheses denote symmetrization of indices. Using Eqs. (21) and (28), we find 
the following gauge transformations for the scalar variables: 

A^A + T' + HT, (29) 

B ^ B-T + L', (30) 

C ^C-d,T + J' , (31) 

S ^ S + HT, (32) 

U ^U + L, (33) 

E^E + d,L + J, (34) 

F ^ F + UT + dJ . (35) 

In a similar manner, vector modes are found to transform as: 



- ■> 



B ^ B + L', (36) 

y -> V + L , (37) 

E^E + d,L. (38) 

In total we have seven scalar modes {A,B,C,S,U,E,F) and three transverse vector 
modes {B,E,V), for a total of ten free functions in Sg^^^. We can now use the four free 
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functions in ^^ to construct six gauge-invariant variables. One possibility is 

^ = A + -[a{B - U')]' , (39) 

^ = S + n{B-U'), (40) 

U = -C + E' + d,iB-2U'), (41) 

A = F + n{B-U')-d,{E-d,U), (42) 
for the scalar sector, and 

f,^K-5,, (43) 

na = Ea-d,Va, (44) 

for the vector sector. Other variables exist, like for example !]„ = £'^ — dzBa- But they are 
not functionally independent; for instance, it is easy to check that E^ = Q'^ + dzTa- 

2. Matter 

We will consider the perturbations to both the infiaton and the anisotropic scalar field. 
Under the gauge transformation (28), any perturbation to -d and ip will transform respectively 

as 

C^^ = Jdz^ , (45) 

C^ifi = ip'T . (46) 

Combining these formulas with the metric potentials, we can define the following gauge- 
invariant perturbations: 

Si} = S^-E + dzU, (47) 

6^ = 6ip + ip'{B - U') . (48) 



Altogether, the system of variables composed by (39-44) and (47-48) span the space of 
cosmological perturbations, implying that fully gauge-invariant linear equations are possible. 
However, in practice it is much easier to perform computations in some specific gauge and 
convert the final result to gauge-invariant expressions. For this purpose, a suitable choice is 

B = E = U = Q = Va. 
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This choice fixes the gauge completely and determines the free functions of Sg^i, and ST^^, in 
terms of gauge-invariant variables. We can therefore perform our calculations in this gauge, 
and translate the final result into gauge-invariant quantities. This is always possible in linear 
theory due to the Stewart-Walker lemma [26]. 

C. Eigenfunctions 

The natural home to define functions on the spatial sections, and in particular scalar 
perturbations, is the Hilbert space of square integrable functions with measure provided by 
the volume element e^dxdydz. Since the Laplacian is a self-adjoint operator on this space, we 
can obtain a spectral decomposition of the identity associated with it. This decomposition 
can then be employed to express any perturbation as a series in terms of eigenmodes of the 
Laplacian. 

In the coordinates defined by Eq. (2) the Laplacian takes the following form: 

v^ = d, + dl + e-^^dl + dl . (49) 

We now want to solve the eigenvalue problem 

V'^(j){x,y,z) = -q'^(j){x,y,z), qeR. 

We recall that the solutions must be (generalized) states in our Hilbert space of square 
integrable functions. The problem can be solved by means of a simple separation of variables, 

(f>{x,y,z) = f{x)g{y)h{z), 

in terms of which the eigenvalue problem becomes: 

1 (df d^f\ -2x^d^g Id^h __ 2 
f \dx dx"^ J 9 dy'^ hdz"^ 

This equation has a nontrivial solution provided that 



ld^g__ ^ ld^h__^ d^f fq'^-n'^ 

g dy'^ ' h dz^ ' du"^ \ u 



2 -^^ ^x^ = -^^ x^ + ^^;;^-i)/ = 0' (^o) 



with real numbers m and n and u = me ^ . The first two equations admit standard plane- 
wave solutions, while the third has a solution in terms of modified Bessel functions: 

f[u) = C^^h{u) + C2V^K,{u) , (51) 
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FIG. 1. Eigenfunctions of the Laplacian in Bianchi III spacetimes. The left panel shows ^/uly 
(continuous line) and y/uK^ (dashed line) for real values of v. The behavior of y/uKu for imaginary 
V is shown in the right panel. The values of {q,n) used were: (0,2), (4,5), and (7,8) (left); (3,2), 
(4,3), and (5,4) (right). 



with the eigenvalue v given by 



V 



1 - 4(g2 



n^ 



(52) 



For real values of z/, the solutions ^/uIy{u) and ^JuKy{u) diverge at m — )■ cxd and u — )■ 
0, respectively. Actually, these solutions are not appropriate for the description of linear 
perturbations, since one can check that they do not belong to the Hilbert space of functions 
under consideration [27] (square integrable with the measure e^dx or, equivalently, du/u'^), 
not even as generalized delta-normalizable states. However, we can find well-behaved (and 
real) solutions if we restrict u to imaginary values. In other words, admissible solutions are 
possible if 

(53) 
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g^ > - + n^ . 
4 

We show some plots of these solutions for some representative values of u in Fig. (1). Note 
that the restriction (53) implies a lower bound on the wavenumber q, which is equivalent 
to an upper bound on the associated wavelength A = 27r/g. If we consider the limiting 
case ra = 0, we see that q > 1/2 = \k\. In other words, at a given instant rj, no physical 
wavelength can be bigger than A7Ta{ri). This is an important signature of the model which, in 
principle, might be detected in CMB multipoles via the Grishchuk-Zel'dovich effect [28, 29]. 
We will however postpone this analysis to a forthcoming paper [30]. 

Returning to the analysis of the eigenfunctions, we have now a complete set in terms of 
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which perturbations can be expanded: 



iu.,m,n)ix,y,z) = ^^^^^^^^e-/'K^i\m\e-)e-ye'-\ (54) 

where u = iu and where the constant prefactor ensures that these functions are delta- 
norniahzed [30]. In conclusion, any perturbation f{x,y,z) can now be expanded as: 

/'OO /'OO /'OO ^^ 

f{x,y,z)= du dn dm (j)^^^^^n){x,y,z)f{uj,m,n) , (55) 

Jo J —OO J —OO 

its inverse being given by: 

f{u, m,n)= / e^'dxdyd^ 0L,m,n)(a;> V, z)f{x, y, z) . (56) 

The symbol * denotes complex conjugation. Note also that, since Ky{u) = K^*{u), only 
positive values of cu contribute to Eq. (55). 

IV. MASTER EQUATIONS 

We are now in an adequate position to identify the effective degrees of freedom and write 
down the master equations for the perturbations. A generic feature of perturbation theory in 
anisotropic spacetimes containing shear is a coupling, already at first order in perturbations, 
between different perturbative modes [18, 31]. Since we are working in a shear-free scenario, 
such couplings will not appear. However, we will still find a residual coupling between 
different scalar modes. In our case, this is due to the way scalar variables are divided in the 
1+2+1 splitting, and also to the spatial curvature of the background spacetime. 

The expressions for the fully linearized Einstein and energy-momentum tensors are pre- 
sented in the Appendix. Here we explain how to gather their combinations into the final 
master equations. 

A. Scalar modes 

All linearized Einstein equations contain scalar modes which can be combined to derive 
master dynamical equations. Before we proceed, note that Eq. (A5) (with a ^ b) implies 
for the scalar modes that 

A = -$ . (57) 
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This extra constraint is similar to what happens in the FRW case when the matter source 
does not contain stress at first perturbative order (which is also the case considered here). 
This relation is always true provided that the spatial average of the perturbations vanishes 
(see Ref. [32]). From now on, we will substitute this constraint in all the equations. 

The first dynamical scalar equation can be obtained from the perturbed Klein-Gordon 
equation. Dropping the tilde in definition (48) for convenience, we find from V^Tq = that 

V' + 2^5¥P' - V^^y? + a^V^^5^ = if' {2^' + AK^ - 2^' - d,U) + 2^"$ . (58) 

Using the constraint (A3) to replace d^H in the expression above, this equation reduces to: 

5ip" + 2H5ip' - (V^ - a'^V^^)Sip = 2ip''^Sip + 2y?"$ . (59) 

The second scalar equation can be found by subtracting (A2) from the trace of (A5). 
We then use the constraint (A3) to substitute S^II in terms of $, \1/, and 6ip. The resulting 
equation is: 

X" + 2nX' - (V^ + Ak)X = AH'^ - 2Lp"bLp . (60) 

Here X = ^ — $. Eqs. (59) and (60) are clearly coupled, and in principle one could 
expect that they were reducible to a single but closed dynamical equation involving some 
combination of ^, $, and ^i^. However, we will show in Subsec. IV C that a closed equation 
for these variables is not generally attainable. Strictly speaking, we can reach an equation 
with those properties only during a slow-roll infiationary phase. 

The dynamics of the scalar sector is also described by the equation for the variable H. 
By combining Eqs. (A2), (A4), (A6), and the corresponding equations for 5T^ we find that: 

H" -t- 2nV^ - (V^ - 2-^' + 4fi;)n = . (61) 

Finally, we have also the dynamics of the perturbed field bd. By combining its Klein-Gordon 
equation with Eq. (A6) one gets (also dropping the tilde in the definition (47)): 

bd" + 2U'd' - V^^ - 4/?^ = . (62) 

It is interesting to note that, even though the scalar sector is being sourced by the two 
matter fields, the above equation is completely decoupled from other perturbed variables 
and, in particular, from the perturbations of the infiaton. This property is reminiscent from 
the fact that ■(? is a massless field without djTiamics at the background level. 
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B. Vector modes 

The dynamics of vector modes are contained in the coupled system of equations (A5) 
and (A6). Because the constraint (A3) does not allow us to decouple them, in principle we 
have to treat Qa and r^ implicitly as one degree of freedom. However, from Eq. (A5) we 
see that: 

dA + f', + 2nTa = f{r],z)Ca, (63) 

where (a is a Killing vector of H^ and / is arbitrary in principle. On the other hand, it is 
not difficult to check that the fact that the perturbations must belong to the space of square 
integrable functions with respect to the volume element defined by the spatial background 
metric implies that the vector (a niust be normalizable with respect to the measure e^dxdy. 
But there exists no Killing vector on EP with this property (see Eqs. (1.4) in Ref. [33]). 
Therefore, the only possibility left is that / = 0. In this case we can use Eq. (63) to decouple 
the system of equations (A5) and (A6), leaving us with one single equation for the variable 

Vt';^ + 2HVl'^ - (V^ + 2K)Vta = Q. (64) 

Remarkably, this equation is formally the same as the equation for the polarizations of 
gravitational waves in a curved FRW universe, where the contribution from the curvature 
appears with opposite sign (see Ref. [24]). This difference in sign is due to the 1+2+1 
splitting of the perturbative modes, as can be checked directly from the perturbed Einstein 
equations (see the Appendix). 

C. Reduced equations 

The system of equations (59-62) and (64) fully characterizes the dynamics. Note that 
Eqs. (61), (62), and (64) are decoupled, and therefore represent three independent degrees 
of freedom of the system. Given that we have eight gauge-invariant variables (6 metric 
potentials plus two scalar field perturbations) and four constraints from Einstein equations, 
we would in principle expect that the coupled scalar equations (59) and (60) could be 
combined into one single scalar equation. When infiation is driven by a single scalar field in a 
fiat FRW background, this reduction is always possible and leads to the dynamical equation 
for the scalar MS variable. In the present case this reduction does not seem possible in 
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its full generality. Curiously, this difficulty does not arise owing to the anisotropy of the 
spacetime, which generically couples perturbative modes, but rather by the presence of the 
spatial curvature. As a matter of fact, precisely the same difficulty appears in curved FRW 
universes, where a nonzero spatial curvature seriously compromises the direct identification 
of the scalar MS variable (see Ref. [34], though). Fortunately, the terms which prevent 
Eqs. (59) and (60) to be combined into one single dynamical equation are precisely the ones 
proportional to the slow-roll parameters. If inflation lasts for at least 60 e-folds, so as to 
solve the standard big bang problems, the dynamics of the scalar perturbations will be given 
by just one scalar MS-like variable. To see that, we introduce a new variable v defined as 



v = a[5'^- ^X 



(65) 



If we now combine Eqs. (59) and (60), using the background equation (12) and the constraint 
(A3), we get the following equation after a lengthy but straightforward calculation: 



v" - Vh - 



z 
z 



2k 






a /^ 
2 \U 



an- 



X 
H 



(66) 



where z = aip' /Ti. 

We now employ the slow-roll coefficients, in terms of which we have 

e'-5'' 



z" 



W 



2 - 3(5 + 2e + (5' - e5 



n 



V2We + 



K e 



5) 



(67) 
(68) 



These expressions are exact. Their expansion at first-order in a slow-roll approximation can 
be directly obtained by means of Eqs. (14) and (15): 

(69) 

(70) 

where in the last equality we have assumed that n/'H? can be neglected compared to e, since 
during inflation l-L = aH grows almost exponentially, whereas n remains constant. Taking 
the above expansions into consideration, we finally obtain: 

1 



1^-1(2 + 66-35) 







7]^ 



-(2 + 6e-3(5) 



0. 



(71) 



This equation is remarkably close to the one found in the isotropic FRW case except, obvi- 
ously, from the different spectrum of the Laplacian. 
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The remaining MS variables are much easier to obtain. With a suitable scaling of 11, -d 
and Qa, Eqs. (61), (62), and (64) can be rewritten as: 



u" - IV' + 2H' + Ak\u = 0, (72) 

e" -(v^- — +4k^9 = 0, (73) 

y:-(v' + "^ + 2K)y, = 0, (74) 
where the new variables that have been introduced are explicitly given by: 

u = all , 9 = aS^ , i/a = aCla . (75) 

Differently from Eq. (71), these equations are valid regardless of the slow-roll approximation. 
For completeness, however, we present their expressions during this stage. Using: 

a" 2 + 3e ^ ,, 1 



r-^ 






U'^-Al + 2e), (76) 



and again neglecting k compared to l-L e ~ e/^ , the above equations become: 



."-^V^ + ^j. = 0, (77) 

^"-(v^-^)^ = 0, (78) 

l/:-(v^ + ^) 1/^ = 0. (79) 

We remind the reader that, despite its appearance, ija is a transverse vector represented by 
one single (scalar) function. 



Equations (71), and (77-79) are the central result of this work, and govern the dynamics 
of the effective degrees of freedom of inflationary perturbations in the considered Bill back- 
ground spacetime. Given that Eq. (71) clearly describes the scalar perturbative mode, and 
that u and ya are pure metric potentials, we can interpret Eqs. (77) and (79) as representing 
the two polarizations of the gravitational waves. We notice that these two polarizations have 
different dynamics, a fact which is a generic feature of anisotropic spacetimes. 

As a final remark, we want to make some comments regarding the massless scalar field 
d used in this work. While its presence is necessary to ensure the FRW dynamics of the 
anisotropic background, its perturbation is completely decoupled from the perturbation of 
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the inflaton field. In other words, equation (71), which is the central equation regarding 
inflationary predictions, remains unchanged regardless of whether we perturb or not i). 
Nonetheless, had we chosen not to perturb this fleld, the equation for u would not be given 
by (72), but rather by: 

u"-(v'-- + 2'hAu = 0, (80) 



a 

as the reader can directly check using the equations in the Appendix and taking ^T'/^s ^ = 0- 
This indirect influence of 6^ in u, as well as the different dynamics of the two polarizations 
of the gravitational waves, might, in principle, lead to observational signatures in their 
spectrum. 

V. CONCLUSION 

The possibility that our current cosmological observations can be made compatible with 
a universe with anisotropic geometry is challenging but really exciting, since this would 
naturally lead us to consider extensions of the cosmological principle beyond our present 
basis for the standard cosmological model. In this context, a especially nice scenario is 
offered by anisotropic spacetimes which are conformally flat (shear- free), since they admit 
an isotropic background expansion, and therefore an isotropic redshift. In this work, we 
have unveiled the dynamics that cosmological perturbations would have in a Bill spacetime 
of this kind during a generic slow-roll inflationary phase. Based on a 1+2+1 splitting of 
the spacetime, we have constructed the full set of linear and gauge-invariant perturbations, 
from which a set of four canonical degrees of freedom have been identifled, and described 
by variables of the MS type in the slow-roll regime. Interestingly, these canonical variables 
have a remarkable similarity to their isotropic (FRW) counterparts. It is worth pointing out 
that, although the identiflcation of the scalar MS variable is only possible during a strict 
slow-roll inflationary phase, the same happens in spatially curved FRW universes. In this 
latter case, the identiflcation of the scalar MS variable becomes severely obstructed by the 
appearance of non-local terms owing to the spatial curvature [34]. We have shown here that 
these terms are of higher order in a slow-roll expansion, and can therefore be ignored during 
inflation. 

In order to expand our equations in natural modes of our system, simplifying in this way 
their resolution and avoiding complicated couplings, we have constructed a complete set 
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of eigenfunctions adapted to the symmetries of the spacetime. This also selects a specific 
functional space that provides the natural home for the perturbations, and allows us to deal 
with issues related with boundary conditions on these perturbations in a straightforward 
manner. A footprint of cosmological perturbations in this cosmological scenario is their 
inability to develop arbitrary super-Hubble modes. This seems to be a general feature of 
spaces with hyperbolic spatial sections, since the same happens in open FRW models as well. 
In this FRW case, analytic extensions of the eigenfunctions are usually invoked in order to 
access the physics of superhorizon perturbations [29] . Whether superhorizon modes will have 
an appreciable impact in future CMB observations or not, and whether such extensions are 
necessary in the present model is a question for further investigation. 
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Appendix A: Useful expressions 

We collect here useful expressions needed to derive the main equations in the paper. Since 
the SV decomposition is unique, different mode contributions can be extracted and treated 
separately in Einstein equations. 

1. Einstein Tensor 

a. Background 



a'G\ = -m'' + l, a'G\=(U^~2-\5l, a^G\ = U^- — + l. (Al) 

a a 



The time index is denoted by 0. 



h. Perturbations 

In the expressions below we keep general the value of the curvature constant k. However, 
one can always replace it with —1/2. 

a^5G% = 6n^^ - 2'HdJl - 2'HA' + D^D^K + Ak^ - AU^' + 2dl^! + D^D^m , (A2) 

o'dGl = D,DyS'^ - ^d^Ta - 2/tf , + ^d^daU - ^d^n'^ + daA' 

-2Hda^ + da<^' , (A3) 

a'^5G\ = -2nd,<^ + 29^^' - ^D^DaU , (A4) 

a^SGl = 51 [2^{H^ + 2H') + 2H^' + dl"^ - ^" - 2H^' - 2Hd,Il - A" - 2nA' 
+D^Dc^ + dl^ - d-Jl' + D^D^A] - D'^^Db)iA + $) + D^^d^VL^) 
+Z^("f ',) + 2'HD^''Th) , (A5) 

a'SG'^ = -9,(9,$ + d,^) + ndaU + ^-d,U' + UdSa + \dS'a 

+^ni' + nn'^ - K^a - ^D^D'^na , (A6) 

a^5G'^ = 2$(?/2 + 2n') + 2n^' + D^D^^ - 2{^" + 2n^' - 2k^) + D^D^^ . (A7) 

As usual, indices between square brackets are anti-symmetrized. 

A useful expression in the derivation of equations for the vector modes is 

DcDaZ = —Za , 
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where Za is any transverse vector. This identity follows from the commutator of two covari- 
ant derivatives and from the fact that Ivj^ = —"jab for the pseudo- sphere. 



2. Energy-Momentum Tensor 

a. Background 

The non-zero components of the anisotropic scalar field ^ and the inflaton field are: 

/^2 /^2 /^2 

'(A9) 



6. Perturbations 

The non-zero components of the perturbed energy-momentum tensor are: 

a^6T^^)o = A-dM, (AlO) 

a2,5T(% , = -n - (50' , (All) 

a25T(':,) , = (A - 9,50) 5,% (A12) 

a25rf^), = 9,50, (A13) 

a^ST^^. = -A + dj(f) . (A14) 



and 



a25T(°^^ = -v?' V + ^v^'^ - a^V'5ip , (A15) 

a^^T"^) , = -ip'da6^ , (A16) 

a'STl>^), = -v'dJv, (A17) 

a2(5T(;) , = (^' V - <l>v^'2 - a V'(5v9)(5," , (A18) 

a^ST^^^ , = V?' V - $v^'^ - a V^v? . (A19) 
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3. Klein-Gordon equations 



a. Perturbations 



The perturbed Klein-Gordon equations for 5d and 5ip can be directly found by means of 
V^T^ = 0. At first order in perturbations, they are given respectively by: 

5d" - V'^6^ + 2H5d' + n' + 2HII - (9^($ + 2^ - A) = , (A20) 

V' - V^^v? + 2H5^' - $>' - 47/^'$ - 2$^" + 2^V' 

+^'d,Il + ^'K' + aH^V^^ = , (A21) 

where V^^p^p stands for d^V/dip^. 
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